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Abstract 
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1. Introduction 



In recent years, the universal scaling in non-equilibrium states have attracted 
much attention [H El IH El- The phase-ordering process (POP) [H] and the 
short-time critical dyanmics (SCD) j2j are two fruitful examples. In the POP, 
the system initially at a very high temperature then is quenched to a heat- 
bath of very low temperature. In the SCD, the heat-bath has the critical 
temperature of the system. 

In both POP and SCD, the initial correlation length is assumed to be 
zero. For finite initial correlation, the scaling invariance of the initial state 
is broken. One would expect a crossover phenomenon that is usually quite 
difficult to study either by theoretic methods or by numerical simulation. 
In order to gain an insight into this phenomenon, in the present paper we 
investigate the kinetic spherical model (KSM) with initial correlation. 

The short-time dynamics of the kinetic spherical model with zero initial 
correlation length has been studied in [HI El EH HB EE21 ECU UH EEH1 EH • In the 
latest publication jUj the effect of non-zero initial order was emphasized. 

The system to be considered in the present paper is initially at an arbi- 
trary equilibrium state. Then it is suddenly put into a heat-bath of another 
temperature. The external field is also removed instantly. The system is 
assumed to evolve following the Langevin equation. We are interested in 
the non-equilibrium transient state in the following time. We will concen- 
trate on the time-dependent order parameter m(t). To include the non-zero 
correlation, one must integrate the non-trivial properties of the initial equi- 
librium state. Though most properties of the spherical model can be found 
in literatures, e.g., the famous book by Baxter [T7j, it is still quite tricky to 
connect the equilibrium initial state with the non-equilibrium dynamics. We 
find that besides the subtraction of mass that is known in the zero initial 
correlation, one need an extra renormalization of the dynamic equation. 

2. The Model 

The Hamiltonian of the spherical model is 

H = f E Sf - (3 J E SiSi - p E (1) 
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with the constraint 

E<s? = * (2 ) 

i 

where < ij > are bonds of a 3- dimensional regular lattice, iV is the total 
number of spins; j3 = -j^f. In the dynamic process, a is a time-dependent 
Lagrange multiplier corresponding to the constraint. 

The fluctuations of spins are defined to be Si = S{— < Si >. In the 
momentum space, one has 

S(p,t) = ^=J2^e tp - ri (3) 



with r\ the position vector of site i. Define 

Q{p) = 1 - ^(cos(pi) + cos{p 2 ) + cos(p 3 )) (4) 

and a function ^ 

w(x) = — -—- — (5) 

The correlation length is inversely proportional to the square root of zo that 
is the solution of the equation 

w(2z ) = (3 J - M ( 6 ) 

for the initial temperature To (/3q = l/ksTo) and the homogenous field h . 
The initial magnetization is given by 

The equation of equilibrium state turns out to be 

, h , 



w( 



' JrriQ 

The initial correlation function is 



&•/(!- m$ (8) 



CM) =< SfeoW.o) >= ^J^ v (9) 
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It is easy to see that 

w(2z ) = ^-Y,C (p, P ) (10) 

p 

The Langevin equation in the momentum space for this constrained spin 
system is 

= _ A(r(t) + (3 jn(p))S(p, t) + r,^ t) (ii) 

where T {(3 = 1/ksT) is the temperature of the heat-bath. The consistency 
condition gives 

r(t) = r su6 + /5j[m 2 (t)-l]+^5:(l-fi(p))<5(-p,t)5(p,t)> 



p 



+ ^E<^(-p.*)»?(p.*)> ( 12 ) 

where the first term comes from the mass subtraction which guarantees 
r c (oo) = at the critical point. Recalling that the equilibrium correlation 
has a zero pole at the critical temperature, one has 

< S(-p,oo)S(p,oo) > \p 



It defines the critical temperature (3 C J = w(0) [TJj. Substituting it into (fT2|) 
at the critical temperature, one can find r su b = 1 for the Ito prescription 
with which the last term of ()12j) is zero due to causality. If the Stratonovich 
prescription is used, the last term of (|T2^) is 1 and = 0. Through the 
paper, Ito prescription will be used. 

By solving (JTTJ), it is not difficult to obtain the response propagator 

i 

i _ -\/3Jfl(p)(t-t')-\ J dt"r(t") 

G P (t, t') = — < S(-p, t) V (p, >= Q(t - t')e *> (13) 

with the Heaviside step function Q(t — t') = 1 for t > t', otherwise Q(t — t r ) = 
0; and the full correlation function (correlation function including the initial 
correlation) 

C p (t,t') = < S(p,t)S(-p,t') > 

= C ( P , P )G p (t,0)G- p (t',0) + C p (t,t') (14) 
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with the correlation function 

oo 

C p {t, t') = 2\J dt"G p (t, t")G p (t', t") (15) 
o 

3. Laplace transformation 

Introducing f(t) = m~ 2 (t) with the time-dependent magnetization m(t) = 
jj < Si >, one can convert the dynamic equation into a linear integrodif- 

i 

ferential equation for f(t) 

^ = 2XPJ - 2\(f3J - l)f(t) + ^ £(1 - n(p))C p (t, t)f(t) (16) 

For convenient, define 7 = (2A/3J) -1 . By Laplace transformation 

00 

F(q) = J dtf(t)e-« 


Equation ()16|) is transformed to 

1 1 7 Co(P,P) 

Substituting Co(p,p) by equation 0, one has 



1 1 7 „ l v 1 

n ' FInJrn 'A N £-< 

F(q) = 



1 - >(7) 

\ + a^2^ Mm) - w ( 2z °)} 



18) 



1 - jl w ^1) 

The last square bracket in the numerator can be written as 

w(r/q) - w(2z ) = [(3 J - w(2z )] - A \pj - w(jq)] (19) 

Where the constant A = 1. In fact, we will see that A plays as a renormal- 
ization multiplier that should be determined self-consistently. The reason 
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will be clear soon. For the infinite system the sum in (jSJ) is replaced by an 
integral. In the continuum-limit one has 

w ( iq ) = W (Q) - D^q) 1 ' 2 (20) 

where the constant D = (^O 3 ^ 2 - In this expansion of w(x), the spatial and 
temporal microscopic detail is lost. However, this microscopic detail would 
have macroscopic effects through equation (fT§j) since it associates with the 
factor 



2z - jq 

which in Laplace reversion is a factor increasing versus time exponentially. 
Therefore, a renormalization multiplier A is needed to rescue the error intro- 
duced by equation (JTUJl . 

Recalling q has the inverse unit of time, one easily finds two characteristic 
time-scales 



. 7 m o 

th '- 



t 



P 



2z 2\ph 
2X^' 2 D " 



where \i = 1 — T/T c . By use of ffT^j) . F(q) is written as 

1 APh 1 



(21) 



F(q) 



yq(l + y/t^q) /3 m% 1 - t h q 

All 

+ ^= 22 

fj,l-t h ql + sft0 



where 



P(l-jjw(2z )) /?J- W (^)l-mg 
When the heat-bath is at the critical temperature, F(q) is 

where B = 



6 



4. Laplace reversion 



Let us first consider the case of T < T c , i.e., \i > 0. 
A direct Laplace reversion to (J22j) gives 



m(t) 



1 + 



+ 



tp — thy th 
A ' h e^erfc 



erfc( 



tp — th 



A 



J/t h \-l/2 



where 



erfc(x) 



dr 



(25) 



(26) 



is the complementary error function. For large x, it has the asymptotic 
expansion 



erfc(x) 



x 



-1/2 



\x-*' 2 + \x~^ + 



(27) 



In equation (|25j) . the last term is not welcome since it leads to a fault exponen- 
tial decay of magnetization. It can be killed by choosing the renormalization 
multiplier as 



A _ 

y/tp ~ \Jth P/i 

One can show that A = 1 in the limit z — > 0. The final result is 

At h 11 



(28) 



m(t) 



1 + 



tp — th\ th 
At h 



erfc( 



erfc( 



{-)Y l/2 



tp — th 

The case of T > T c can be attained by similar method 



(29) 



m(t) 



At h 



t 



tp — th V th 



th 



+ 



1 + 



At t 



tp — th 



^{2-erfcU-))}- 1/2 



(30) 



One can attain the critical evolution of magnetization from the Laplace 
reversion of (|24|) . 



IB 2 



1 -V4 r 




1-1/2 



m(t) 



(31) 



7T 



where 



poo 



e 



-y 



g(x) = dy 



(1 + xy) 1 / 2 



(32) 



.70 



It is clear that g(0) = 1. The critical behavior also can be attained by taking 
the limitation T — > T c in equation (|2T?|) . 

To recover the results of zero initial correlation of [Oj, one only need to 
take the limitation of (3q — > oo, that is zo — > oo (or ^ — > 0), in (j2S| , (JSOJ) and 
(13 1|) . respectively. In this limit, A = with defined in jHj and t/j = Ttt^. 

5. Discussions and Conclusions 

In summary, we have studied the transient behavior of KSM that is quenched 
from an arbitrary temperature into another. The formula can describe or- 
dering/disordering phenomena and the critical dynamics. We find that the 
correct long-time behavior can be only recovered after the subtraction of 
mass as well as the renormalization of the dynamic equation. 

From the exact magnetization obtained in the present paper, one can 
find out the flows of renormalizational group of the bare parameters j3, zq 
and mo (or f3o) under the scale transformation. Changing the time scale by a 
factor b z , in order to keep the macroscopic quantity m(t) unchange up to a 
scaling factor, th and tp must be transformed in the same way, 



The relations of bare parameters and the scaling factor b are implicitly defined 
in the above equations. These are the so-called characteristic functions for 
the dynamic crossover phenomena |H1 El UH UHj • 




(33) 



while A must be an invariant, 



A(b) = A(l) 



(34) 
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